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EFFECTS OF CIRCULAR GEOMETRY IN SIMULATION OF CONVECTION 


IN ROTATING-SPACECRAFT TANKS 


E. Dale Martin, Barrett S. Baldwin, and Margaret Covert 
Ames Research Center 


SUMMARY 


For computational simulation of the convection and mixing of stratified fluid in a rotating 
tank (such as used in Apollo flights) with time-dependent rotation, the Navier-Stokes convection 
problem is formulated for a circular-tank configuration. The final equations result from a general 
approximate theory derived previously for combined forced and contained natural convection in a 
time-dependent rotating system. The equations are cast in terms of vorticity and stream function in 
a form most convenient for computation, with a convenient transformed coordinate system, and 
appropriate boundary conditions are derived for use in subsequent numerical computations. 
Accurate representations for the cryogenic supercritical oxygen thermodynamic properties are used 
in the computations, and a highly efficient numerical finite-difference scheme and computational 
method, used previously by Lomax and Bailey, are employed. 

Computations in both square and circular two-dimensional tank configurations simulate the 
motion and heating of the oxygen in spherical tanks under conditions corresponding to Apollo 
flights. The circular-tank calculations are performed to compare with previous square-tank results 
and to confirm the previous general conclusions. The effects of rotation reversal maneuvers, spin-up 
maneuvers, and heater cycles on the evolution of temperature distributions and on the pressure and 
potential for pressure collapse are examined and compared for the two configurations. 

The results show that appreciable mixing can be achieved in a short time by small rotation 
maneuvers to prevent pressure collapse. The conclusions from previous calculations of a square-tank 
configuration are confirmed. 


INTRODUCTION 


This report describes extensions of a previous study (ref. 1 ) of combined forced and natural 
convection and mixing of a thermally stratified fluid in a spacecraft tank with time-dependent 
rotation. 

The background of the study is described in some detail in reference 2, but for completeness 
here a description of the physical problem is provided below. The study applies to spacecraft 
fluid-storage systems in general (including, e.g., those in orbiting space stations); but the study arose 



in connection with the oxygen tanks in the Apollo spacecraft, so the particular conditions imposed 
and the specific results obtained refer to Apollo oxygen tanks. 

In spacecraft, oxygen and hydrogen are stored in a supercritical cryogenic state because a large 
mass of fluid can be stored at high density (small volume) and at pressures low enough for reason- 
ably lightweight tanks and high enough to ensure adequate supply to various systems. In Apollo, the 
oxygen tanks are spheres with inside diameters of 63.5 cm (25 in.) and distances from their centers 
to the spacecraft axis of 91.44 cm (3 ft) and 152.4 cm (5 ft). (The 3-ft rotation arm is used in the 
computations to be described later.) Initially, the tanks are loaded with 150 kg (330 lb) of liquid 
oxygen. Then the oxygen is heated by an internal electric heater to vaporize the oxygen. The heater 
is then operated periodically to maintain a system design pressure of 620 X10 4 ± 24X10 4 N/m 2 
(~900 ± 35 psia). When the pressure drops to the lower design limit because of withdrawal of fluid, 
the heater is automatically turned on. When the pressure rises to the upper limit, the heater is 
switched off. These heater cycles and pressure cycles are discussed in later sections of the report. 

Under these conditions, the oxygen temperature can be of order 150° K, near the critical 
point. The fluid is supercritical and is not distinguishable as either a liquid or a gas. At these 
conditions, a small change in temperature can cause a substantial pressure drop. If the fluid is 
significantly stratified in density and temperature and is suddenly mixed, such a resulting pressure 
drop (also called pressure collapse or pressure decay) may even drive the fluid into the undesirable 
two-phase state and may render inoperative the systems supplied by the oxygen-storage system. For 
understanding of these concepts, the reader should refer to Kamat and Abraham (ref. 3). 

During portions of flights by spacecraft such as Apollo, the conditions can develop for which 
pressure collapse can occur. When the craft is in a zero-gravity condition with no accelerations, each 
heater cycle results in increasingly severe temperature gradients near the heater. The gradients decay 
very slowly because of the inefficiency of pure heat condition for energy transfer. Therefore, to 
prevent the potentially dangerous levels of stratification, frequent mixing must be provided. Prior to 
the Apollo 13 accident, the fluid inside the tanks was circulated by means of electric fans. The 
Apollo 1 3 Review Board found that the electric fans played a part in the accident that aborted that 
mission, and the oxygen tanks for subsequent missions were redesigned to eliminate the fans. An 
analysis determined that the stratification problem could be crucial, so other means of stirring, or 
inducing convection, were needed. One problem was to determine the adequacy of vehicle rotation 
and maneuvers to produce sufficient convection and mixing through the body forces due to 
temperature differences and rotational accelerations. 

As a means for estimating the mixing effectiveness of rotation and rotational maneuvers, the 
fluid convection can be simulated by numerical computations. For this purpose a set of equations 
was developed (ref. 4) that represent a general approximate theory for the Navier-Stokes description 
of fluid convection in a rotating system, including the influence of both the effective buoyancy 
body forces due to temperature and density stratifications and the arbitrary time-dependent rota- 
tion of the tank. The equations reduce to a generalization of the Boussinesq approximation in the 
special case of steady rotation with variable buoyancy body forces, but also properly represent the 
forced convection due to time-dependent rotation. The equations are valid for combined forced and 
contained natural convection, whereas previous treatments of convection in rotating systems had 
dealt with either pure natural convection or pure forced convection (e.g., see ref. 5, sec. F, 22; and 
refs. 6 and 7). 
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For convenience in the use of the numerical methods employed (and in keeping with the spirit 
of LighthilPs description of the flow development and its relation to vorticity ref. 8, pp. 57—60), 
the equations were cast in terms of the vorticity and stream functions in reference 4. For the first 
numerical calculations (refs. 9 and 10), convection in a rotating square two-dimensional tank was 
treated. Although three-dimensional effects such as axial flow due to Ekman-layer suction (e.g., 
ref. 6) and development of Taylor vortices (see ref. 7) could not be represented in a two- 
dimensional formulation, the first numerical treatment by this theory (refs. 9 and 10) yielded 
significant information regarding the adequacy of vehicle rotation and of rotation reversal for 
mixing. 

The first numerical treatment of the Apollo oxygen tank convection employed the square-tank 
configuration for convenience in developing the numerical methods, since many of the essential 
two-dimensional effects could be simulated. A two-dimensional circular tank, with a small circular 
interior boundary at the center and a heat source near the center, may more closely simulate the 
actual spherical oxygen tank containing the central circular-cylindrical capacitance gage and the 
heater, as described in reference 2. Therefore, one purpose of this report is to extend the previous 
study to circular geometry — i.e., to formulate, in a form particularly suited to subsequent numeri- 
cal computation, the problem of convection in a circular tank and to compare the results of 
computations for the circular geometry with those for the square geometry. The flow equations 
that are used result from the general approximate theory (based on a small-density-variation 
approximation) developed in reference 4 for combined forced and contained natural convection, 
with influence by both buoyancy body forces and time-dependent rotation of the vehicle. The final 
equations for computation are put into forms very similar to those used for the square tank in 
reference 4. -The accurate representations of thermodynamic properties from reference 1 1 are used, 
as well as the special procedures from reference 12 for analysis of the thermodynamic states 
resulting from the small-density-variation approximation. 


FLOW EQUATIONS AND PROBLEM DEFINITION FOR CONVECTION 
IN A ROTATING CIRCULAR TANK 


Configuration and Flow Equations 

The circular-tank configuration for analysis of the convection is shown in sketch (a). The r, 6 
polar-coordinate system is fixed relative to the tank, with 6 measured from a line perpendicular to 
the rotation arm (from the tank center to the center of rotation, of length R c ). The radius vector 
from the rotation center O to a point r, 6 in the tank is denoted by R, with components R r and Rq 
in the directions of the unit vectors ,e r and (the directions of increasing r and 6 , respectively): 

R= e r R r + e d Rq (la) 

where 

R r = r + R c sin 6 (lb) 

Rq = R c cos d (1 c) 
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Sketch (a). Circular tank configuration. 


shown in sketch (a), a, may be prescribed by its 
a, so that 


The time-dependent angular velocity £2 is 
positive when counterclockwise. Then in the 
r,d coordinates the effective body force per 
unit mass, given by equation (2.12) in refer- 
ence 4 is 


£ 

= £ r 2r + 

&e s$ 



(2a) 

where 





Sr 

= R r n 2 

+ 2f2 vq 

+ £IRq 

~a r 

(2b) 

se 

= R d n 2 

- 2 £2v r 

- nR r - 

~ a d 

(2c) 

in 

which ar 

and aQ 

are the- components 

of a 


prescribed vector a,, which is an arbitrary 
time-dependent linear acceleration of point O 
in the noninertial system of coordinates fixed 
on the tank, relative to an inertial system. (As 
an effective body force per unit mass, £ may 
be taken to include any body forces acting on 
the fluid in the inertial reference frame.) As 
magnitude | aj and instantaneous relative direction 


& ~ &r a r + &d a 6 

(3a) 

a r - !aJcos(a -6) 

(3 b) 

a.Q = |a,|sin(a - 0) 

(3c) 


Also in equations (2), £1 is defined to be dSl/dt, and v r and v$ are the respective r and d com- 
ponents of the fluid velocity vector V measured relative to the coordinate system fixed on the 
rotating tank: 


V = &r v r + &d v d ( 4 ) 

The approximate system of flow equations developed in reference 4 is based on the full 
compressible Navier-Stokes equations with the assumption that the thermodynamic state variables, 
density p, pressure p, and temperature T, have small variations in space and relatively slow 
variations in time. The approximation proceeds by a small-perturbation approach referenced to the 


4 



thermodynamic state ( p 0 , p 0 , T 0 ) representing nearly constant values of p, p, and T that are 
averages over the space of the tank at a given time. For the small-perturbation analysis, one write's 


P - P 0 (* + 5 ) 

where 

5 * -0(7’- T q ) |8| «| 


Along with the above assumptions, the specific enthalpy is taken to be independent of pressure; 
and the transport properties — viscosity p = p Q v and thermal conductivity k — are assumed to be 
nearly constant in space and slowly varying in time. That is, spatial gradients of p and k are 
neglected, but p, k, and the specific heat c p are all allowed to vary with time and to depend on the 
spatial distribution of temperature and on the average pressure at any time. 

To satisfy the approximate equation of conservation of mass resulting from assuming 
|5|«1 (V V = 0), one can define a stream function (eqs. (2.33e) in ref. 4) such that 

l dip d\p 

Vr = ~r He ’ v o = ~ Yr (6) 

It is advantangeous to use as a primary variable the vorticity (Lighthill, reference 8, pp. 57-60), 
which is defined by 


(5a) 

(5b) 

(5c) 


w=v x v 

where for two-dimensional flow (with the unit vector = ji r X / e^): 

co = e 7 oj 

and the scalar co is (cf. ref. 4, eqs. (2.63) and (2.64)): 

> 

13 1 9>T 



(7a) 


(7b) 


( 8 ) 
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From reference 4, eq. (2.58c), the approximate equation for the rate of change of vorticity, derived 
from conservation of momentum, is 



bco 1 1” 9 

a 

^ r 9 / 

r 9o>\ 

co> ”1 



IF * 7 [aF ( "V“ ) + 

ae (|, 8“) 

= a- 2 i S' \ 


+ 90 2 J 





a t 

/ 1 

9 T\ 





set Yr + 


To)-™ 

(9) 

The approximate energy equation, with Dp/Dt neglected (ref. 4, eq 

|. (2.84) 

with (2.85) and (2. 75)), is 

bT 1 

9 9 

/ k \ 

.if 9 / 

9 T\ 

9 2 r"| v 


02 1 

+ 

^ 1 


5- 

n 

1 ^ r 97 V 


90 2 J + c p $ 

(10a) 


where the dissipation function <f> is 


= 2 


r / av A 2 ( * .vi r 9 ( v q \ 1 9 .i 

L \ br ) + \ r 90 + r ) J + |_ r br ^7 J + 7 90 J 


(10b) 


Application of Boundary Conditions 


Initially (at t = 0) a temperature distribution is specified according to a prescribed temperature 
stratification; and v r -vQ= 0 everywhere (assumed-rigid-body rotation) so that 

co = 0 everywhere initially. (1 1 ) 


Boundary conditions for t > 0 must be specified on \p, cu, and T for application to equations (8), 
(9), and (10). 

On the boundaries r - a (outer tank wall) and r - (interior circular boundary representing 
the capacitance gage in the Apollo oxygen tank), the conditions of no normal mass flux and no slip 
are 


or 



(12a) 



(12b) 
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where the subscript b indicates a value at a boundary. The application of these conditions to the 
Poisson equation (8) for \p and the parabolic equation (9) for to must be formulated in a manner 
to be compatible with the numerical-solution technique. Such application is described as follows. 

For an arbitrary smooth two-dimensional boundary with coordinates s and n measuring 
respectively distances along and perpendicular to the boundary, and with curvature K = K(s) of the 
boundary surface, the conditions are 


-1 3 \}j 

no mass flux: = - — — — — = 0 on the boundary (13a) 

n 1 + Kn 3s 

3i p 

no slip: v„ = — = 0 on the boundary (13b) 

a on 


where and v n are the respective velocity components tangent and normal to the boundary. Again 
with the subscript b denoting a boundary value, equation (13a) can be replaced by 


no mass flux: ~ c(t) 

and equation (13b) is 


no slip: 



04) 


(15) 


Also in terms of the coordinates s and n, it can be shown that conditions (13) along with V 2 'p = -cu 
give 



(16a) 


At a small distance e measured perpendicular to a smooth boundary, Taylor’s series for \p = \p e is 
so that insertion of equations (15) and (16a) gives the approximate condition 


to / 


-20 Pe ~ 'Pb) 


(16b) 


Equations (14) and either (15) or (16) can be used to represent (13). 
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In most problems, \pb can be specified. In that case, the conventional approach is to specify 
\pb =0 (eq. (14)) for application to equation (8) for the <// field. Effectively, condition (13a) has then 
been applied, and \p is known everywhere at the current time step of the computation. Then 
condition (16b) is specified as a boundary condition at the current time step in solving the vorticity 
equation by an explicit scheme for interior values of cj at the next time step. Effectively the no-slip 
condition, (13b) or (15), has then been applied through (16b) to the vorticity equation. 

It may happen, as is the case in the present problem at the interior boundary, that \pb is not 
known. One could apply (14) to equation (8) and iterate on the value of c{t) at each time step, but 
this would consume machine-computing time. Therefore, a suggested alternative is to apply the 
no-slip condition (15) to the Poisson equation (8) for i jj\ then when i p is known at the current time 
step, effectively use the no-mass-flux condition incorporated into the Taylor series that results in 
condition (16b) (with \pb constant at a given time) for application to solving the vorticity equa- 
tion (9). (This alternative approach is believed to be new.) 

From the above discussion, for the circular-tank problem one can now specify the following 
boundary conditions on \p and co for application to equations (8) and (9), respectively, 




r - a : =0 




-2(^-Ar 

(Ar) 2 


r = a. 


d\}/ 

dr 


= 0 


co =co fli = = 


+Ar-*a i y 

(A rf 


(17a) 

(17b) 

(18a) 

(18b) 


The application of (17) and (18) is illustrated in sketch (b). 



„ No slip, ^ - 0, applied to V 2 tf/ = - w 
or 

\ „ No mass flux, included in co - <o a(I (eq. (18b)), 

applied to vorticity equation 

' - - No mass flux, = 0, applied lo V 2 vJ/ = - u 

— No slip, included in co-w a (eq. (17b)) , 
applied to vorticity equotion 


For application to equation (10a), the 
heat flux or dT/dr can be specified on all 
boundaries (with value zero for insulated 
boundaries). A heat source at some interior 
point near r -a x can be specified to represent 
the power input from the heater. 


In d, all variables must be periodic with 
Sketch (b). Application of boundary conditions. period 27T. One variation in the numerical 

application of the boundary condition on 
vorticity is discussed on page 13 of the sec- 
tion on Computational Procedure and in the Discussion section, page 1 8. The above conditions 
could be applied without alteration if a sufficiently fine computation mesh were used. 
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The function £2(t) must be prescribed and, in the example of sudden rotation reversal, can be 
represented by 

•' 1 «(*■) = « 0 t-l + 25(/)| •- (19a) 

Cl(t) = m o S(t) (19b) 

where is a constant angular velocity, S(t ) is the unit step function, 

5(0 = 0 t < 0 (20a) 

= 1 t> 0 (20b) 

and S(t) is the Dirac delta function. Simple numerical approximations to equations (20) may be 
used, such as a linear change of 5(f) from 0 to 1 between t =0 and small time t l , for which the 
corresponding approximation to 5(0 is the constant value \/t j for 0 < t <0 and zero outside that 
time interval (used also in ref. 4). 

Transformation of Variable and Dimensionless 
Equations for Computation 

To facilitate numerical computations using a “Poisson solver” for the computation step involv- 
ing solution of equation (8) for i//, it is convenient to make the transformation 

r = ln^ (21) 

which puts the Laplacian operator into a form similar to the rectangular-coordinate form. 

Then with L as an arbitrary length, and A T as an arbitrary constant reference temperature 
difference, the following dimensionless variables and parameters are introduced: 
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£ = £(?) 


= j=(j)e^ r = — 
L \L / L 2 


c cL 2 T - T q \ p 

co*= — , H= , * = - 

v A T v 


VqL 


VyL 

u = — , v = 

v v 


2fif * 2n o 

EL *(r)= , ft =- — 

vL~ 2 ° vL ~ 2 


Gj = Gr 


Sr 




Gi — Gr 


£0 


a o 2L 


a,* = , a 2 

n 0 U 


a 0 




(22a) 


<j>* = — (j> 

p 2 


with: 


Reynolds number = Re 


<V 2 

v 

P 0 vc p 

Prandtl number = Pr - — ; — 


Grashof number = Gr - 


Eckert number = Ec = 


( EL 0 2 L)L 2 ^T 

V 2 

(EL n L) 2 


c p AT 


(22b) 


(This particular nondimensionalization is not motivated by any order of magnitude considerations 
for the present problem but was used for convenience in comparing test calculations with results 
obtained from similar previous calculations.) Equations (6), (8), (9), and (10) become, in dimen- 
sionless form (dropping the stars from ca*, EL*, <f>*, a, *, and a 2 * from here on): 



1 


H 

£ 


( a// 

-G 2 g^ + G, 


a h\ 

30 ) 


dEl 

dr 


(23) 
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M + _L 1 

3t £ 2 _3£ 


_a_ 

30 


J = _L /J_\ (pH + &Ji\ Ec_ 
J * U 2 / \ 9 ^ 2 M 2 / + Re 2 


( 24 ) 


a 2 ^ a 2 * 

ar 2 dd 2 5 


(25) 


f/ = I ^ F= zli* 

£ 30 ’ \ ar 


(26) 


where 


Gi - — — (| + sin0 ) t ^ 2 2 + fi F + (cos0 ) ^ 

Re 2 4 2 ar 


Gr 

fle 2 


[■ 

[ 


(cos0) jft 2 -«(/-(£ + sin0) ^ ^ 
4 2 (2T 


•-[(f *(H’]*H 


3F 

a? 


-Gr a x 

(27a) 

- Gr a 2 

(27b) 

acA 2 

dd I 

(28) 


and where £ is the function of f defined in equations (22a). (In eq. (24), the time derivative of T 0 in 
H has been neglected.) 


The initial conditions for the dimensionless variables in equations (23)through (26) are 


r = 0: U=V= co = 0 

H (f,0, 0) prescribed 


(29) 


The transformation of variables can be 
thought of, and the equations analyzed, in 
terms of a rectangular f , 6 mesh as illustrated 
in sketch (c), corresponding to the r,d polar 
coordinates shown. 

The boundary conditions in the £0 
plane include the periodic condition for all 
dependent variables: at each f in sketch (c), 
all dependent variables must have the same 
values at 6 = 2ir as,at 0 = 0. 



Sketch (c). Transformation of coordinates. 
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The boundary conditions at the outer tank wall corresponding to conditions (17) are 


f = 0 : tf=0 


'2(*>r = o - iAr i 


CO = 


(a/L) 2 (AO 2 

At the interior circular boundary, conditions (18) become 

a . 9 ^ 

ir° 


(30a) 

(30b) 


Cl) = 


+iArr ( ^ = ? fll ] 


(ajL) 2 (Af) 2 


(31a) 


(31b) 


for application to equations (25) and(23) respectively. For specified heat flux - k(dT/dr ) at r - a 
and r = a, , the corresponding values of 



(32a) 


(32b) 


are specified for application to equation (24). 

The dimensionless function H(r) and d£l/dr may be specified, for example, as in the discussion 
of equations (19) and (20), with t and t x replaced by r and Tj . 


COMPUTATIONAL PROCEDURE 


The numerical method for the integration of equations (23) through (25) is very similar to that 
described in reference 9, which should be referred to for details. An explicit finite difference 
procedure is used for equations (23) and (24) that employs the predictor-corrector method of 
reference 13. The main difference in the overall procedure is that the stream-function equation in 
reference 9 was solved at each step in the integration over time by a Buneman recursive-cyclic- 
reduction technique, whereas equation (25) with the periodic boundary condition at 6 = 0 and 
6 = 2n and with conditions (30a) and (31a) is solved in the same context in the present method by a 
Fast-Fourier-Transform technique as described by Hockney in reference 14. The thermodynamic 
relations and attendant computational procedures are described in references 1 1 and 12. Boundary 
conditions on the dimensionless relative temperature H are imposed corresponding to a specified 
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heat leak through the walls. A time-dependent value of H is imposed on computational grid points 
selected to represent the heater used to overcome pressure decreases due to removal of fluid from 
the tank. 

Although pressure gradients must be considered in the momentum equations, pressure 
gradients are negligible in the energy equation for the low-velocity cryogenic oxygen flows under 
consideration. In the energy equation the variation of pressure with time has only a cumulative 
effect that does not affect the motion of the fluid. As discussed in reference 4 the only important 
coupling between the energy and momentum equations is through temperature gradients and con- 
vection of temperature variations. As a result, for the purpose of finding the motion of the fluid, 
the density can be considered a function of temperature with pressure as a slowly varying time- 
dependent parameter that is uniform in the tank. Time-dependent temperature and density distribu- 
tions then result from integration of the equations of motion over short time steps. Computation of 
the evolution of pressure and other thermodynamic quantities is described in reference 1 2. One such 
quantity of particular interest is the potential pressure decay, which is the drop in pressure that 
would result from complete adiabatic mixing of the stratified fluid. Thus, the potential pressure 
decay is a measure of the pressure loss that could result from spacecraft maneuvering after a given 
level of stratification has built up from nonuniform heating of the fluid in a near zero-g environ- 
ment (see, e.g., ref. 3). 

Several expedient remedies for removal of difficulties arising from application of the numerical 
method are discussed in references 10 and 12. A uniform spacing of computational grid points was 
employed. It was found that it is not possible to resolve adequately the viscous boundary layer 
that develops near the wall. One could use the boundary conditions (17b) and (18b) precisely if a 
sufficiently fine mesh were used. However, to avoid excessive computing time a coarse mesh (17 X 17 
computation points) was employed in most of the calculations. It then was found necessary to put 
limits on the absolute value of vorticity so that extreme (unreal) viscous stresses would not be 
propagated to the interior from the walls when large flow velocities occurred near the walls. Thus, 
the boundary conditions (17b) and (18b) were overridden and replaced by constant values when the 
computed values exceeded the fixed upper limit. Most of the calculations used an upper limit equal 
to four times the vorticity corresponding to the solid body fluid motion associated with a spacecraft 
rotation reversal with a total change of 1 .6667X10" 3 revolution per sec (6 revolutions per hour). 
Sample calculations were made with an upper limit half as large to establish that variations of 
thermodynamic quantities, such as pressure and potential pressure decay, were insensitive to this 
approximation. The effect of the approximation on the results is considered further in the 
discussion section of this report. 


RESULTS 

Rotation Reversal Maneuver 

Calculations were made for both square and circular tanks in which the spacecraft was taken to 
be rotating initially at 8.333X10" 4 revolutions per sec (3 revolutions per hour). An assumed initial 
stratified state was imposed in each case with temperature varying linearly across the tank. The 
warmer fluid was placed in the stable position toward the center of rotation. As shown in 
reference 10, this linear type of temperature distribution does occur after prolonged heater cycling 
during rotation at a steady rate. In the absence of spacecraft maneuvering other than constant 
rotation, the calculations showed a very slow decrease in potential pressure decay as expected. The 
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decrease in this case is due to heat conduction arising from the mild temperature gradient and 
negligible motion of the fluid with respect to the tank. When the direction of rotation of the 
spacecraft was reversed, however, a swirling motion of the fluid ensued as shown in figure 1 . 

Photographs are presented of a cathode-ray 
display tube on which the velocity vectors are 
plotted at a number of points in the field. In 
the square tank, the velocity vectors are 
shown at the computational grid points. In 
the circular tank, the vectors are plotted at 
points in a rectangular array that does not 
correspond to computational grid points. This 
was done to avoid nonlinear concentrations of 





Figure 1 . Velocity vectors for square and circular tanks 
showing clockwise swirl introduced by rotation reversal. 


vectors in the display that would be difficult 
to view. Linear interpolation was used to 
obtain the values at the display points. 

The upper and lower photographs on the 
left in figure 1 correspond to a time 360 sec 
after the rotation reversal. The photographs 
on the right show the velocity vectors 480 sec 
after reversal. The off-center swirls visible in 
these frames move continuously in a 
clockwise direction around the center in the 
same direction as the general motion of the 
fluid. The flow patterns in the square and 
circular tanks are similar as can also be seen in 
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Figure 2. Temperature distributions after rotation 
reversal (square tank). 


figures 2 and 3. The flow velocities are of 
order 0.003 m/sec (~0.01 ft/sec) near the outer 
boundaries. 



720 sec 1320 sec 

Figure 3. Temperature distributions after rotation 
reversal (circular tank). 
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Figure 2 shows a series of temperature distributions in a square tank. For temperatures greater 
than the average, the deviations from the mean are represented by vertical lines. As an aid to 
visualization, for temperatures less than average the deviations from the mean are shown as hori- 
zontal lines. In either case, the length of the line indicates the magnitude of deviation from the 
mean. The upper left photograph in figure 2 shows the initial assumed temperature distribution, 
which varies linearly from top to bottom. The warmer fluid (vertical lines) is toward the center of 
rotation, which is about 1-1/2 tank widths (diameters) below the center of the tank. The remaining 
photographs in figure 2 show the temperature distribution at later times after the rotation reversal. 
In the lower right photograph considerable mixing and dissipation of the temperature distribution 
have taken place. Figure 3 shows a similar evolution of temperature distributions in a circular tank 
resulting from a rotation reversal maneuver. 

Figure 4 shows a comparison of the 
decreases in potential pressure decay (defined 
on p. 13) resulting from rotation reversal of 
square and circular tanks. These calculations 
were based on the van der Waals approxima- 
tion of the thermodynamic properties of 
oxygen described in reference 12. No large 
difference in the degree of agreement 
exhibited for the two geometries is to be 
expected with the use of more accurate 
thermodynamic relations. No heater opera- 
tion or heat leak were allowed for in these 
calculations and the pressure dropped to 
about 4.8X10 6 N/m 2 (~700 psi) near the 
ends of the runs. For such adiabatic mixing 
the decrease in potential pressure decay is 
equal to the drop in actual pressure. The initial linear temperature distribution for the circle was 
chosen to cover a larger range than that for the square so that the initial potential pressure decays 
would be equal. These results show that no large differences in the response of the fluid state to 
rotation reversal maneuvers are to be expected for square or circular geometry. 



Figure 4. Comparison of reductions of potential pressure 
decay resulting from rotation reversal in 
square and circular tanks. 


Heater Cycles 

As discussed in references 10 and 12 and elsewhere, the Apollo oxygen tanks are equipped 
with heaters to overcome the drop in pressure from removal of fluid. Pressure sensors in the tanks 
turn the heaters on when the pressure falls below about 6.0X10 6 N/m 2 (~870 psia) and turn them 
off when the pressure rises above about 6.3X10 6 N/m 2 (~910psia). During fluid removal, the 
heaters are thus automatically cyclically turned on and off about every 1 200 sec, depending on the 
flow rate and the amount of fluid remaining. Figure 5 shows typical velocity and temperature fields 
that develop after several heater cycles for different simulated tank geometries and heater positions. 
The flow field in the upper left photograph was generated with elements on the left wall acting as a 
heater. A steady spacecraft rotation rate of 1.1111X1 CT 4 revolutions per sec (0.4 rph) was imposed 
in all of the cases shown in figure 5. In the upper right photograph the velocity vectors resulting 
from an internal position of the heater are shown. In the lower left photograph, the velocity vectors 
generated in a circular tank with an internal heater position are shown for comparison. The lower 
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Velocity Temperature 


Figure 5. Velocity and temperature distributions resulting 
from heater operation with uniform rotation of 
the spacecraft at 1.1111X1 CT 4 rps (0.4 rph) 
(hexagon denotes position of internal heater). 


right frame shows the temperature distri- 
bution in the tank corresponding to the same 
calculation. In times of order 7200 sec (2 hr) 
with small spacecraft rotation rates, only a 
small volume of fluid near the heater is heated 
appreciably. It can be seen in the lower right 
photograph of figure 5 that the heated fluid 
tends to move toward the center of rotation, 
which is about 1-1/2 tank diameters below 
the center of the tank. 

Figures 6, 7, and 8 show the pressure 
cycles and potential pressure decay versus 
time corresponding to the geometries of the 
upper left, upper right, and lower photo- 
graphs in figure 5. In these calculations the 
accurate thermodynamic relations of Stewart 
(ref. 15) were used as described in refer- 
ences 1 1 and 12. In figure 6, the pressure rises 
when the heater is on and subsides when the 
heater is turned off. The lower plot of 
figure 6 shows that the potential pressure 
decay responds similarly, but also undergoes a 
cumulative rise as the amount of stratification 
increases. Comparison of figures 6, 7, and 8 
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Figure 6. Pressure cycles for square tank with wall 
heater, rotation rate 1.1111X1 CT 4 rps 
(0.4 rph). 


Figure 7. Pressure cycles for square tank with internal 
heater, rotation rate 1.1 1 11X1 0 -4 rps 
(0.4 rph). 
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shows that the response of pressure and 
potential pressure decay to heater cycling is 
insensitive to the geometry and position of 
the heater in the tank. The notation 
“Quantity 90%” in the figures indicates that 
the. tanks contained 90 percent of the initial 
mass in these calculations. “Heat Leak 10 W” 
refers to the rate of heat leakage into the tank 
through the walls. 


IO^N/m^ psia 




Figure 8. Pressure cycles for circular tank,rotation rate 
1.1111X1CT 4 rps(0.4rph). 


l0 6 N/m* psio 


Spin-Up Maneuver 

Figure 9 shows the effect of vehicle 
spin-up on the pressure and potential pressure 
decay. In this calculation the spacecraft was 
not- rotating initially and the cumulative rise 
of potential pressure decay was greater than 
in the previous figures, where the rate of 
rotation was 1.1111X1 Or 4 rev/sec (0.4 rph). 

The spin-up was started at 7620 sec after the 
start of the calculation and was completed in 
about 300 sec. The rotation rate was held con- 
stant at 8.333X1CT 4 rev/sec (3 rph) thereafter. 

After a delay of 600 sec the swirling motion 
of the fluid introduced by the vehicle spin-up 
had a large effect on the pressure cycle. The 
delay can be attributed to the fact that the 
initial response of the fluid was essentially a 
solid body rotation relative to the tank. A 
period of time elapsed before the interaction 
of density gradients with the new acceleration 
field resulted in nonuniform motion of the 
fluid. The resulting mixing action caused a 
drop in pressure and also reduced the poten- 
tial pressure decay. Thus, even though the 
heater was turned on at 8940 sec, when the 
pressure fell below 6.00X 10 6 N/m 2 
(870 psia), it continued to drop because of the 
mixing action and was down to 
5.58X10 6 N/m 2 (810 psia) at the end of the 

calculation. The calculation ended at that time because of failure to allow for lower pressures in the 
machine code as described in reference 1 2. A similar drop in potential pressure decay and pressure 
resulting from spin-up after heater cycling was obtained for a square tank in reference 10, showing 
that the effect is insensitive to the geometry (square or circular) used in the calculation. These 
results also show that the potential pressure decay associated with local hot spots produced by 
heater cycling is more effectively removed by spacecraft maneuvering than a linear temperature 
distribution across the entire tank. 



Figure 9. Effect of spin-up on pressure cycle (circular 
tank, initial rotation rate 0, final rotation 
rate 8.333X1 O' 4 rps or 3 rph). 
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This can be seen by comparison of figures 4 and 9. In figure 9 the onset of an appreciable 
decrease in potential pressure decay occurs more quickly and the decrease is more rapid. 


DISCUSSION 


The foregoing results show that appreciable mixing can be achieved in a short time by changing 
the spacecraft rotation rate by modest amounts. No significant differences were found in the 
calculations for circular geometry in this report and those for a square tank in reference 10. As 
discussed in reference 10, to obtain numerical stability in the flow near the corners of a square tank 
it was found expedient to put an upper limit on vorticity. The upper limit imposed was four times 
the value corresponding to the solid body rotation of the fluid associated with a rotation reversal 
from 8.333X1CT 4 rev/sec (3 rph) to 3 rph in the opposite direction. This limit on the absolute value 
of vorticity also resulted in slip flow in the boundary layer near the walls. It was argued that the 
resulting obscuration of possible reversed flow eddies in the corners is not objectionable since the 
actual tanks being simulated are spherical and do not have corners. It was also demonstrated by 
estimates based on flow over a flat plate that the boundary layer at the walls would remain 
relatively thin during times simulated by the computations. Therefore, the slip flow at the walls 
resulting from the imposition of an upper limit on vorticity would affect the motion of only a small 
fraction of the total volume of fluid. It was concluded that the calculations for a square tank 
provide conservative estimates of the mixing induced in the main bulk of fluid by the interaction of 
density gradients with the acceleration field due to spacecraft maneuvering. Nevertheless, it was 
deemed worthwhile to make calculations for a circular tank to determine whether the methods used 
to suppress corner flow effects could cause errors that had not been recognized. (In the circular- 
tank calculations, the limiting values of vorticity did not occur except right on the boundaries.) 

The present calculations for circular tanks in which no large differences were found relative to 
the earlier calculations for square tanks reinforces the conclusions derived from the square-tank 
calculations. The levels of potential pressure decay to be anticipated according to our calculations 
are in reasonable agreement with previous estimates from Apollo 1 2 data and stratification analyses 
(see ref. 1 0). We concur in the conclusion that dangerous levels of stratification will not occur with 
mass flow rates and acceleration fields corresponding to past and anticipated Apollo operational 
practices. 


Ames Research Center 

National Aeronautics and Space Administration 
Moffett Field, Calif. 94035, June 28, 1972 
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